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Abstract 


Using coarse-grained molecular dynamics simulations, we analyze me¬ 
chanically induced dissociation and unfolding of the protein complex 
CD48-2B4. This heterodimer is an indispensable component of the im¬ 
munological system: 2B4 is a receptor on natural killer cells whereas 
CD48 is expressed on surfaces of various immune cells. So far, its 
mechanostability has not been assessed either experimentally or the¬ 
oretically. We find that the dissociation processes strongly depend on 
the direction of pulling and may take place in several pathways. In¬ 
terestingly, the CD48-2B4 interface can be divided into three distinct 
patches that act as units when resisting the pulling forces. At experi¬ 
mentally accessible pulling speeds, the characteristic mechanostability 
forces are in the range between 100 and 200 pN, depending on the 
pulling direction. These characteristic forces need not be associated 
with tensile forces involved in the act of separation of the complex be¬ 
cause prior shear-involving unraveling within individual proteins may 
give rise to a higher force peak. 
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1 Introduction 


A variety of processes in biological cells are performed and controlled by large 
protein complexes rather than single proteins. Examples include replication, 
transcription, translation, signaling, protein transport and degradation, and 
cell adhesion [T] . Contemporary structural biology and molecular biophysics 
are thus concerned with understanding of ever larger and more dynamically 
complex molecular assemblies. 

Force spectroscopy is a powerful method to study mechanical properties of 
single biopolymers and biomolecular assemblies. A typical force spectroscopy 
experiment uses atomic force microscope (AFM), or optical tweezers, to mea¬ 
sure the response of a biomolecular system to stretching, squeezing, or tor¬ 
sional manipulation M0U51IM- These experiments have provided a vari¬ 
ety of information on transport and force generation in cells, DNA replication 
and transcription m, DNA unknotting and unwinding |Sjj, ATP synthesis 
in mitochondria [9], mechanical behavior of virus capsids under squeezing 
mm, and other biologically relevant problems. 

There is a growing recent interest in theoretical and experimental assess¬ 
ments of mechanical stability of protein complexes. If there are no covalent 
bonds between individual proteins then stretching by the AFM tip may re¬ 
sult in dissociation of the complex, unraveling of the constituent proteins, 
or both. The outcome reflects the strength of the binding forces and the 
structural organization of the complex. Examples of the systems considered 
so far include the titin ZlZ2-telethonin m nu, dimeric tubulin ca, the 
3D domain-swapped cystatin PHHH, oligomeric fibrinogen [IS], covalently 
and non-covalently linked homo-dimeric systems with the cystine knot motif 
E2B El], and the homo-trimeric capsomers of the CCMV virus capsid ra- 
In particular, Dima and Joshi na have identified a viscoelastic behavior in 
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unraveling of the intermeric couplings of tubulin chains. In distinction to the 
fibrous systems [TB], there is no physiologically motivated direction of pulling 
in the cystine knot, capsomer, and cystatine complexes [201 EH E2UI6UIZ] - 
In these systems, the response to stretching has been found to be strongly 
anisotropic: It depends on which site the tip is attached to and which site 
serves as an anchor. This is a significant extension of the anisotropy phe¬ 
nomena that have been observed in monomeric proteins [T9l |. 

In this paper, we direct our attention to the anisotropic effects arising 
in mechanical manipulation of heterodinrers. As an illustration we consider 
a cell adhesion complex in which mechanical manipulation may lead either 
to unfolding, or to dissociation, or to both effects combined, depending on 
the geometry of pulling. Specifically, we focus on the murine CD48-2B4 
complex [23] with the Protein Data Bank structure code of 2PTT. Protein 
2B4 is a heterophilic receptor present at the surface of natural killer cells. 
The antigenic CD48 proteins are expressed on surfaces of various immune 
cells. Interactions of CD48 with 2B4 regulate the activity of the natural 
killer cells, which are components of the innate immunity against tumors 
and virally infected cells. Even though the CD48-2B4 complex appears not 
to have been studied through the single molecule manipulation techniques 
yet, it serves as a convenient system for theoretical epxlorations. 

The crystal structure of the CD48-2B4 complex that we use in this study 
comprises two extracellular immunoglobulin-like domains: the C2-type do¬ 
main of CD48 (shown in blue in Fig. [TJ), and the CD48-binding domain of 2B4 
(shown in red in Fig. Q]). These domains comprise 110 and 112 amino acid 
residues, respectively, which means that the CD48-2B4 system is too large 
for a thorough all-atom molecular dynamics study on relevant time scales. 
Here, we use a coarse-grained structure-based model to study the response 
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of the CD48-2B4 complex to stretching at constant speed. For simplicity, 
we consider stretching only by the terminal sites. Some ways of choosing 
them may result in partial or complete unfolding of the constituent proteins, 
and other may lead to splitting of the complex into separate parts. We ex¬ 
plore the roles of the inter- and intra-protein interactions in the competing 
processes of unfolding and dissociation of the CD48-2B4 complex. 

The 2B4 receptor binds to CD48 with a modest strength, with the equi¬ 
librium dissociation constant in the low micromolar range and very fast on 
and off rates j22j. In the native state, 19 residues of 2B4 are in contact with 
16 residues of CD48, and most of the contacting residues belong to loops and 
turns in CD48 and in 2B4. The CD48-2B4 interface involves only 36 con¬ 
tacts, as determined by atomic overlaps (see the Methods section), and one 
of our focal points is quantitative analysis of how and when these contacts 
break during the dissociation process. 

We find that the dissociation scenarios as well as the force-displacement 
patterns depend sensitively on the choice of the termini that are used to 
implement stretching. In some cases, the proteins separate without any no¬ 
ticeable deformations or structural changes; in other cases, CD48 unfolds 
partially before the proteins dissociate. However, the inter-protein contacts 
always break in groups rather than sequentially. The 36 contacts between 
CD48 and 2B4 can be divided into three distinct patches, which are denoted 
here as a, b, and c, as shown in Fig. |J] We find that the three interface 
patches cooperate as units in resisting the stretching forces. 

We also find that a significant resistance to stretching arises due to two 
types of mechanical clamps, which involve shear and tensile forces, respec¬ 
tively. In shearing, two /3-strands slide by, see the left-hand panel of Fig. [21 In 
tensile strain, the separation between two residues in contact increases along 
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the line connecting them, as illustrated in the right-hand panel of Fig. [21 
Some choices of the pulling termini generate merely tensile effects, followed 
by separation of the protein complex. Other choices involve both types of 
mechanical clamps, in which case several force peaks arise on the path to 
separation. The largest of the force peaks needs not be tensile in origin, 
and its height is not necessarily related to the binding strength of the com¬ 
plex. This observation is relevant for the interpretation of force spectroscopy 
experiments on protein complexes. 

Since there is no force spectroscopy data on the CD48-2B4 complex that 
we could compare with our simulation results, we also simulate stretching of 
synaptotagmin 1, which is a membrane-trafficking multidomain protein that 
has been studied thoroughly in single-molecule pulling experiments [21 [22 
125] . Synaptotagmin 1 facilitates fusion of synaptic vesicles. Its cytoplasmic 
part is composed of two domains, C2A and C2B, which together form the 
C2AB complex. Our simulation results agree with experimental findings. In 
particular, we find that C2AB is significantly less stable than the 127 domain 
of titin, and that the C2A domain in the complex is less stable mechanically 
than its C2B partner. 

2 Methods 

2.1 Coarse-grained model 

We use our coarse-grained structure-based model mmmm, in which 
each amino acid residue is represented by a single bead centered on its C Q 
position. The beads are tethered together into chains by strong harmonic 
potentials with the spring constant febond = 100 e/A 2 , where e is the depth 
of the potential well associated with the native contacts, which serves as the 
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basic energy scale in our model. The native contacts are identified using 
an overlap criterion [30] applied to the coordinates of all heavy atoms in 
the native structure. However, the amino acid pairs that are very close 
sequentially, (i,i + 1) and (i,i + 2), are excluded from the contact map [28] . 
It should be noted than our scheme to identify native contacts is different 
than in, for example, self-organized polymer model [15l [18] where a uniform 
length cutoff is used. 

The interactions within the native contacts are described by the Lennard- 
Jones potential 


H NAT (r^) = 4e 




( 1 ) 


Here, ry,- is the distance between residues i and j in the contact, and the 
parameters are chosen so that each contact in the native structure is sta¬ 
bilized at the minimum of the Lennard-Jones potential. The value of e is 
approximately given by (110 ± 30) pN A, as has been determined by com¬ 
paring simulational results to the experimental ones on a set of 38 proteins 
[29] . The contacts between the proteins are treated in the same manner as 
the contacts within the proteins as both sets are dominated by hydrogen 
bonds. The interactions in the non-native contacts are purely repulsive and 
given by the truncated and shifted Lennard-Jones potential corresponding 
to &ij = ro/ \/2 with ro = 4 A. The energy function comprises also harmonic 
terms that favor the native values of local chiralities in each amino acid chain 
[•TTI . A harmonic potential with the same force constant fcbond as the C Q -C a 
pseudo-bonds is used to model the disulfide bond between Cys 3 and Cys 100 
in the 2B4 receptor (chain B in the PDB structure 2PPT). 

The solvent is implicit and the system evolves in time according to the 
Langevin dynamics. The overall force acting on a particular bead i is a sum 
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of three terms: (i) the direct force that derives from all the potential terms, 

(ii) the damping force that is proportional to the velocity of the bead, and 

(iii) the random force, T,., that represents thermal noise. The corresponding 
equations of motion 


m 


d 2 r; 
d t 2 
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( 2 ) 


are solved by the fifth order predictor-corrector algorithm with the time step 
of 0.005 r. Here, 7 is the damping coefficient and all beads are assumed to 
have the same mass m. The dispersion of the noise is given by i/2y ksT , 
where fc# is the Boltzmann constant and T denotes the temperature. All 
simulations were performed at ksT = 0.3 e which is near-optimal in folding 
kinetics and is of order of the room temperature. The damping coefficient is 
set to 7 = 2m/ t. This value corresponds to the overdamped case - practically 
Brownian dynamics - and the characteristic time scale, r, is of order 1 ns, 
as argued in Refs. [321 E3|. 

The native contacts can break during the time evolution. Our criterion 
for this to happen is that the distance between residues i and j in a contact 
exceeds 1.5 07 ,-. 


2.2 Pulling simulations 

Stretching of the CD48-2B4 complex is implemented by attaching harmonic 
springs to two terminal amino acids. One of the springs is fixed in space 
and the other one is moved with a constant speed, v p , so that the distance it 
travels in time t is d = v p t. The force constant of the pulling springs is taken 
as K = 0.12e/A 2 , which corresponds to about 1 pN/nm and is close to the 
elasticity of typical AFM cantilevers [22J. 

In our simulations, the response force, F, acting on the pulling spring 
is measured and averaged over time periods that correspond to the spring 




displacements of 0.5 A [55], The F-d curves (see the top panels in Fig. [3] and 
and several later figures) may come with several peaks, and the height of the 
largest of them will be denoted by F max . We perform simulations for a range 
of speeds with the trajectories evolving for up to 4 ms. Specific examples of 
the F-d patterns will be shown for v p = 0.002 A/r. 

All of the pulling simulations start from the native state. In the course of 
the simulations, the breaking and re-formation of native contacts is followed 
in time. The native contact between residues i and j is considered broken if 
the inter-residue distance rjj exceeds a cutoff length, as defined in section IjTTl 
Due to thermal fluctuations, the broken contacts may get re-established. To 
characterize the unfolding and dissociation pattern of the events, we record 
the spring displacements at which the native contacts break for the last time. 
In the corresponding scenario diagrams (see the bottom panels in Fig. [3] and 
several later figures), native contacts are labeled by the sequential distance 
\j—i\ between the residues involved. Here, we use the convention in which the 
index i runs from 1 to 216 so that the values between 1 and 106 correspond 
to CD48 and between 107 and 216 - to 2B4. Positions of five amino acids 
are not available in the structure hie (residues 1, 2 and 3 in chain A, and 
residues 1 and 110 in chain B) and hence the effective sequential size of the 
system is smaller than the nominal 106 plus 110. The inter-protein contacts in 
patches a, b, and c (shown in Fig. Q]in orange, purple and green, respectively) 
correspond to i around 30 and j around 145, i around 40 and j around 195, 
and finally to i around 90 and j around 150, respectively. 

Our model is based on the knowledge of the native structure of the pro¬ 
tein complex. Both CD48 and 2B4 are /3-proteins, see Fig. [D Protein 2B4 
comprises two apposing /3-sheets. One of them is formed by five /3-strands 
(/3i, /3s, fa, /3 3 and ^ 4 ) and the other one by three /3-strands (/3 2 , /3q and 
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(3$). Here, the /3-strands are labeled from the N-terminus to the C-terminus. 
2B4 also contains two short helices: a\ located between /3 4 and fa, and a 2 
between fa and fa. There are also two /3-sheets in CD48: the first one is 
formed by fa, fa, fa, fa, fa and fa, whereas the second one is composed of 
fa, fa and fa. Protein CD48 contains also a short a-helix, denoted here by 
cci, which is located between fa and fa in the apposing /3-sheets. 

3 Results and Discussion 

For any heterodimer, there are six distinct ways to choose a pair of termini, 
and each choice may lead to more than one response pathway in stretching. 
Here and below, the termini in CD48 will be denoted by N and C whereas 
those in 2B4 by N' and C', see Fig. [Q 

Stretching is applied along the line connecting the two chosen termini. 
To describe the shift in the orientation of stretching in the protocols dis¬ 
cussed below, we introduce an angle 6 between the C-C' direction and other 
directions of stretching. We obtain 6 ) (N-N , )=68°, 6 I (N-C , )=16°, 6{ C-N')=23°, 
0(N-C)=28°, and 0(N , -C')=38°. 


3.1 Dissociation 

We first discuss situations in which the two termini defining the pulling di¬ 
rection belong to different protein chains. This way of pulling results in 
mechanically induced dissociation of the two proteins. 

3.1.1 N-N' protocol 

It is useful to first consider the N-N' protocol, in which the pulling springs 
are attached to the N-termini of CD48 and 2B4. We observe existence of two 
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kinds of dissociation pathways, as illustrated in Fig J3] In both of them, there 
is no breaking of contacts within 2B4 because the disulfide bond between 
Cys 3 and Cys 100 holds the chain in a sturdy manner. 

On the first pathway, protein CD48 unfolds partially through unzipping 
giving rise to a minor force peak at the spring displacement d of about 55 A. 
This force peak is associated with breaking of contacts between the 
and the fa-ai loops, and between the (3 1-/?2 loop and the C-terminal tail of 
CD48. The unzipping process stops at d of about 100 A, when strands /3i 
and /?2 in CD48 get fully unraveled. At d between about 100 A and 120 A, 
the protein complex re-orients itself so that a tensile clamp between the two 
proteins is formed, see the right-hand panel of Fig. [2j Then, at d of about 
125 A, all of the interface contacts break almost simultaneously, which yields 
large response forces with F max of 2.2 e/A, which is comparable to that of the 
127 domain of titin [29], 

On the second pathway, there is no appreciable unfolding within neither 
of the two units. Instead, the two force peaks that are observed are due to 
rupturing of the interface contacts. The contacts in patch b rupture at d of 
about 55 A, which is reflected in moderate response forces with F max of about 
1.3 e/A. Then the contacts in patches a and c rupture almost simultaneously 
at d of about 70 A with a smaller force. In both rupturing events, we observe 
tensile mechanical clamps. What prevents unraveling of CD48 on the second 
pathway are the contacts between j3\ and /3g] these contacts break early 
on the first pathway but they never do on the second pathway. The first 
pathway occures more often but the level of dominance varies with the speed 
of pulling, as shown in the top panel of Fig. SI in Supporting Information 
(SI). 
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3.1.2 N-C' protocol 


Pulling in the N-C' protocol also results in two pathways that lead to sepa¬ 
ration of the protein complex, see Fig. [4j The two pathways are very similar 
initially, at d < 80 A, when and /3 2 unravel. However, in the first pathway, 
F max of about 1.6 e/A arising at d ~ 80 A is due to shear between /3 2 and 
/3g, and between /3 2 and (3 7 in CD48, as indicated by arrows in Fig. [4] The 
inter-protein patches break later: first patches a and c at d ~ 100 A and 
then patch b at d of about 120 A. The associated force peaks are smaller. 

In the second pathway, unraveling of fd\ and /3 2 gives rise only to a minor 
force peak at d ~ 80 A. The next peak at d of about 110 A is due to 
rupturing of contacts in patch a. However, F max of about 1.8 e/A arises 
when the contacts in patch b and those between /3 3 and /?§ in CD48 break 
simultaneously. These events occur at d of about 130 A and are immediately 
followed by the breakage of contacts in patch c at d ~ 140 A. The first 
pathway dominates and its statistical weight increases with decreasing speed 
v p , see the middle panel of Fig. SI in SI. In fact, the second pathway is never 
observed at small pulling speeds ( v v < 0.001 A/r). 

3.1.3 C-N' and C-C' protocols 

In protocols C-N' and C-C' we observe only tensile effects and no multiple 
dissociation pathways. The corresponding F-d patterns and dissociation sce¬ 
narios are shown in Fig. [5j For both protocols, no intra-protein contacts are 
ruptured before the two units break loose. In the C-N' protocol, the contacts 
in patch b get broken at d of about 30 A, which yields F max of about 1.3 e/A, 
and then the contacts in patches a and c break simultaneously at d of about 
45 A with smaller forces. In the C-C' protocol, the three patches break 
simultaneously at d of about 35 A, and F max increases to about 1.9 e/A. 
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3.2 Unfolding: N-C and N'-C' protocols 

We now consider situations in which the two termini defining the pulling 
direction belong to the same protein chain. Pulling in the N'-C' protocol 
results in a monotonic growth of F with d because of the tethering effect 
exerted by the disulfide bond between Cys 3 and Cys 100 in 2B4, see the 
dashed line in Fig. [5j On the other hand, pulling in the N-C protocol results in 
full unraveling of CD48 that can proceed along two pathways. The maximum 
force arises either as the first (pathway 1) or the third (pathway 2) peak in 
the F-d curves, as indicated by arrows in the top panels in Fig. |6j The 
values of F max are 1.7 e/A and 1.9 e/A, respectively. In the first pathway, 
the mechanical clamp is due to shearing between (3\ and /? 9 in CD48, as 
illustrated in the left-hand panel of Fig. [21 The shearing between /3i and /3g 
occures at d of about 30 A. In the second pathway, the mechanical clamp 
is due to shearing between f3$ and /3 9 , and between /3 3 and (3%, which takes 
place at d & 140 A. 

Interestingly, the larger the sequential distance \i — j\ of a contact, the 
earlier its rupture takes place. The contacts in the inter-protein patches 
get broken only after the tertiary structure in CD48 is destroyed. Each 
of the patches acts coherently, but unlike what happens in the dissociative 
protocols, individual patches get affected at different times. 

We observe the same two pathways of CD48 unraveling in the presence 
and absence of the partner protein 2B4, as indicated by the solid and dashed 
lines in Fig. [6J respectively. Interestingly, as shown in the bottom panel of 
Fig. SI, the first pathway is statistically more relevant at intermediate pulling 
speeds (0.0002 A/r < v p < 0.005 A/r), and the second pathway dominates 
at low and high pulling speeds. 


13 


3.3 Shear and tensile mechanical clamps 

As discussed before, F max arises either as a result of shear between /3-strands 
in CD48 or due to the tensile clamps between CD48 and 2B4, depending on 
the pulling protocol and pathway. The distinction between the two kinds of 
mechanical clamps sits in the directionality of the movements of the residues 
involved, as illustrated by the simulation snapshots in Fig. [2j However, it 
also shows in the time dependence of the distances between the residues in 
contacts that form the clamp. The left-hand panel of Fig. [7] shows three 
such distances as a function of d in pathway 1 in the N-C protocol. These 
distances correspond to three native contacts between f3\ and /3g in CD48. 
For d < 32 A, the inter-residue distances r i3 fluctuate around their native 
values. At a critical distance of about 32 A, shearing between the two /3- 
strands leads to the simultaneous breaking of the three contacts, which is 
reflected in the fast increase in the inter-residue distances r t] from about 
5 A to about 35 A. In the right-hand panel of Fig. [TJ the three distances 
rij correspond to contacts between CD48 and 2B4 in patches a, b and c, 
respectively. The situation considered here relates to pathway 1 in the N-N' 
protocol. For d < 125 A, the inter-residue distances r l3 fluctuate around their 
native values. At a critical distance d ~ 125 A, the contacts in patches a, b 
and c break simultaneously, which leads to a sudden increase in the distances 
from about 10 A to about 70 A. The difference between the two clamps is 
in the steepness and magnitude of the increase of the inter-residue distances 
- from 5 to 35 A, and from 10 to 70 A in the case of shear and tensile clamps, 
respectively. The simulation snapshots in Fig. [2] correspond, left-to-right, to 
the panels of Fig. |T| 
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3.4 The velocity dependence 


Our discussion so far has been focused mainly on trajectories obtained at 
one particular speed of pulling, v p = 0.002 A/r. We now discuss how F max 
varies with v p , see Fig. HI We performed simulations at seven different speeds 
of pulling ranging from 10 _4 A/r to 10~ 2 A/r, which corresponds to the ex¬ 
perimental speeds between 10 nm/rns and 1 nm//xs. At each of the pulling 
speeds we ran at least 60 trajectories. We separated the simulation data 
into pathways, if any, and averaged over the trajectories. We find an efective 
logarithmic dependence |34j ES] 


T max = P In (v p /vo) + q (3) 

in the range of the speeds studied. Here, the unit of speed is u 0 =l A/r. The 
fitting parameters p and q are summarized in Table Q] The parameter p de¬ 
scribes how fast F max changes with v p . The range of variations in p is seen to 
be within a factor of 1.5, which reflects anisotropy effects. Interestingly, the 
two pathways of the N-C protocol exhibit almost identical values of the pa¬ 
rameter p. The slowest dependence is observed for the C-N' protocol whereas 
the most rapid one is seen for the second pathway of the N-N' protocol. 

Assuming the Bell-Evans model [STj [35j in which 

k B T , (x^Ku 


F = 

1 rna.Y 


ad \k B Tk 0 


(4) 


one can extract the intrinsic off-rates of the interacting proteins, ko, and 
the location of the free-energy barriers confining the proteins in their bound 
state, ad, from the fitting parameters p and q. Comparison of Eqs. (ED and 
gD gives ad = k B T/p and k 0 = e~ q ^ p Kvo/p. The resulting values of ko and 
summarized in Table |U In case of deviations from the logarithmic 


x* are 


dependence gD, one can use more sophisticated theories [36j that permit to 
extract also the height of the free-energy barriers AGA 
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It is interesting to observe that the smallest off-rates in Table U correspond 
to either complete (protocol N-C) or partial (N-N' pathway 1) unfolding of 
CD48. The off-rates that correspond merely to dissociation (protocols C-C', 
C-N', and N-N' pathway 2) are orders of magnitude larger. Nevertheless, 
some values of k 0 as given in Table [I] seem unreasonably small. They have 
been obtained at T = 0.3 e/fcg. Higher temperatures should lead to their 
substantial increase. Note that the room temperature is more likely to corre¬ 
spond to about 0.35 e/ks (notice that the value of e itself comes with about 
25 % uncertainty [25] ) and the physiological temperatures are still higher. 

We note that when the pulling sites are the N- and C'-termini, the second 
pathway occurs only at the pulling speeds v p > 0.001 A/r, which is indicated 
by the absence of data points at small speeds v p < 0.001 A/r. To further 
investigate this observation, we calculated the frequency of occurrence of 
the first pathway in protocols N-N', N-C' and N-C. At a particular speed of 
pulling, the frequency zq is defined as the ratio of the number of trajectories 
at which the first pathway has occured to the total number of trajectories. 
The occurrence frequency zq as a function of the pulling speed v p is shown in 
Fig. SI. In the N-N' protocol, the frequency zq varies between 0.6 and 0.95, 
and attains the maximum value at v p = 10~ 3 A/r. In the N-C' protocol, zq 
increases from 0.65 at v p = 10 -2 A/r to 1 at v p — 10~ 3 A/r, and remains 
constant, zq = 1, in the interval between v p = 10~ 3 A/r and v p = 10 -4 A/r. 
In the N-C protocol, the frequency zq increases with the speed v p in the 
interval between v p = 10 -4 A/r and v p = 10 -3 A/r, and decreases with v p in 
the interval from v p = 10~ 3 A/r to v p = 10 -2 A/r, with the maximal value 
of 0.7 at v p = 10” 3 A/r. The error bars in Fig. SI indicate the standard 
deviation of the Bernoulli process. 

We now turn our attention to the process of CD48-2B4 separation. The 
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separation event takes place at the moment when no contacts are present 
between CD48 and 2B4 for the first time on the trajectory. The correspond¬ 
ing spring displacement at which separation takes place will be denoted by 
d s . Fig. S2 in SI shows that d s decreases with v p in a way that depends on 
the protocol and the pathway. Its most significant dependence is observed 
for pathway 2 in protocal N-C', which is the pathway that finds no continu¬ 
ation at small speeds. The decrease in d s reflects the greater role of thermal 
fluctuations at lower speeds since fluctuations foster separation. 

The spring displacement d does not directly relate to the end-to-end dis¬ 
tance, L, between the pulling termini. In fact, L reflects transformations in 
the structure in a more direct manner. Fig. S3 in SI shows that L s - i.e. the 
value of L corresponding to separation - also depends on v p , although usually 
less sensitively. For instance, there is much less variation in L s compared to 
d s in the C-C' protocol. On the other hand, they are about the same for 
both pathways in the N-C' protocol. 

3.5 Stretching of the multidomain synaptotagmin 1 

To our knowledge, there have been no single-molecule pulling experiments 
performed on the CD48-2B4 complex. Thus the assumptions underlying 
the construction of our model, especially about the relative strength of the 
inter-protein contacts remain untested. To establish connections to force 
spectroscopy experiments on membrane-associated multidomain proteins, we 
performed pulling simulations of synaptotagmin 1. The latter protein is 
involved in remodeling the plasma membrane during neurotransmitter release 
at the synapse. The cytoplasmic region of synaptotagmin 1 contains two 
domains, C2A and C2B, which together form the C2AB module. These 
two domains are quite similar structurally but have been found to exhibit 
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different mechanical stabilities within the C2AB module [21, 07]. 

The issue which is analogous to that pertaining to the CD48-2B4 complex 
is whether it is indeed sensible to assume that the inter-domain contact 
energy, which we denote here by e', is about the same as the contact energy e 
within each domain. Fig. E] shows a typical F-d curve and the corresponding 
scenario diagram for e' — e and v p = 0.002 A/r. The highest force peak 
(labelled as B in Fig. [9]) is associated with rupturing the core of the C2B 
domain. The second highest peak (labeled as A in Fig. E]) arises when (3i 
and @2 of the C2A domain unravel. This event is followed by breaking the 
C2AB interface, which results in a minor force peak (labeled by I in Fig. [9]). 
A typical simulation configuration at this stage is depicted in the lower right 
corner of Fig. [HI It is similar to the relevant configuration obtained in steered 
molecular dynamics simulations performed at a pulling speed that is three 
orders of magnitude larger than used here 

After extrapolation to the experimental pulling speed, see Fig. El we 
obtain the values of F max to be equal to (100 ± 30) pN and (70 ± 20) pN for 
the force peaks B and A, respectively. These values are consistent with the 
experimental results of about 100 pN and 50 pN (with the statistical error 
of about 20 pN) for the C2B and C2A domains [24] , 

We now consider the situation in which e' is distinct from e. We find that 
if e' = 2e then the C2A domain becomes much more stable than the C2B 
domain, which disagrees with experimental findings. On the other hand, if 
e' = e/2 then the second highest force peak becomes associated with ruptur¬ 
ing the core of the C2A domain rather than with unraveling of f3\ and /? 2 , 
see Fig. S4 in SI. However, the characteristic forces (1.6 e/A and 1.75 e/A 
for C2A and C2B, respectively, at v p = 0.002 A/r) are very close to those 
presented in Fig. [9] (F^ x = 1.6 e/A and F ® ax = 1.8 e/A) at the same speed 
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of pulling. Nevertheless, the appearance of the whole F-d pattern is much 
closer to the experimental one when e' — e than when e' is reduced. We 
therefore conclude that assuming the equality of e 1 and e is consistent with 
the experimental results and with the similarity of the physical nature of the 
interdomain contacts. 


4 Conclusions 

The mechanical stability of adhesion proteins and their complexes is neces¬ 
sary to sustain interactions and signaling between cells. To investigate the 
mechanical stability of the adhesion protein complex CD48-2B4, we have per¬ 
formed pulling simulations using our coarse-grained structure-based model. 
We find that the force peaks arise either as a result of shear between /9-strands 
in CD48 or due to tensile mechanical clamps between CD48 and 2B4, de¬ 
pending on the mode of pulling and dissociation pathway. If there are several 
force peaks in a given F(d) curve then, in general, the maximal one may be 
associated with shear or with the tensile separation. In the former case, the 
measurement of F raax would not provide an estimate of the adhesion forces. 

Table fTOl provides a summary of our findings regarding the values of F max 
for various protocols of stretching. For a given v p , the table is symmetric. We 
present it in a way which gives values of F max for two pulling speeds: above 
the diagonal the data correspond to v p of 5 x 10 -3 A/r - the standard used in 
our surveys of mechanostability ra - and below the diagonal they correspond 
to 6 x 1 CF 6 A/r, which is in the lower range of speeds used in protein- 
related AFM experiments m- The highlighted diagonal blocks in the table 
correspond to same-chain pulling which results either in structure unfolding 
(N-C) or just deformation (N'-C'). The off-diagonal blocks correspond to 
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pulling of different chains which results in dissociation. The v p dependencies 
are seen to be comparable in these two classes of pullings, and dissociation 
may come either with a smaller or larger F max than unraveling. 

We observe existence of strong anisotropies. For dissociation, F max varies 
between 0.9 and 1.7 e/A at the extrapolated experimental speed. For unrav¬ 
eling - between 1.35 and 1.55 e/A. In the N-C protocol, and also on the first 
pathway of the N-C' protocol, F max is caused by shear between /3-strands in 
CD48. In other pulling protocols, F max is due to tensile mechanical clamps 
between CD48 and 2B4. It is important to note that the largest force signal 
during dissociation may actually come from shearing and not tensile separa¬ 
tion. 

Both CD48 and 2B4 are anchored in membranes through the C-terminal 
ends [23]. The forces acting on the CD48-2B4 complex during cell adhesion 
are thus applied in the C-C' direction (although membrane fluctuations on 
nanometer scales may alter this direction to some extent [38] )• Interestingly, 
we find that in the C-C' pulling simulations neither CD48 nor 2B4 unravels 
(both proteins maintain their structure) and the entire CD48-2B4 interface 
acts as a whole to resist external forces with F m ax = 1.4 e/A at the experi¬ 
mental speed of pulling. The latter value corresponds to about 150 pN. These 
results suggest that the CD48-2B4 complex is adapted to resist significant 
forces while maintaining the structure of individual proteins. In particular, 
the disulfide bond between Cys 3 and Cys 100 prevents the rupturing of 
contacts within 2B4. 


20 


Acknowledgments 

This work has been supported by the Polish National Science Centre Grants 
No. 2012/05/B/NZ1/00631 (BR) and 2011/01/B/ST3/02190 (MC) as well 
as by the ERA-NET Grant FiberFuel (MC). LM acknowledges access to the 
KRUK computer cluster at the Faculty of Physics, University of Warsaw. 


References 

[1] Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. Molecular 
Biology of the Cell. New York: Garland Science; 2002. 

[2] Carrion-Vazquez M, Oberhauser AF, Fisher TE, Marszalek PE, Li H, 
Fernandez JM. Mechanical design of proteins studied by single-molecule 
force spectroscopy and protein engineering. Prog Biophys Mol Biol 
2000;74:63-91. 

[3] Rittort F. Single molecules experiments in biological physics: methods 
and applications. J Phys Cond Mat 2006;18:531-583. 

[4] Vogel V. Mechanotransduction involving multimodular proteins: con¬ 
verting force into biochemical signals. Annu Rev Biophys Biomol Struct 
2006;35:459-488. 

[5] Nauman KC, Nagy A. Single-molecule force spectroscopy: optical 
tweezers, magnetic tweezers and atomic force microscopy. Nat Meth 
2008;5:491-505. 

[6] Crampton N, Brockwell DJ. Unravelling the design principles for single 
protein mechanical strength. Curr Opin Struct Biol 2010;20:508-517. 


21 


[7] Galera-Prat A, Gomez-Sicilia A, Oberhauser AF, Cieplak M, Carrion- 
Vazquez M. Understanding biology by stretching proteins: recent 
progress. Curr Opin Struct Biol 2010;20:63-69. 

[8] Witz G, Stasiak A. DNA supercoiling and its role in DNA decatenation 
and unknotting. Nucl Acids Res 2010;38:2119-2133. 

[9] Buzhynskyy N, Sens P, Prima V, Sturgis JN, Scheuring S. Rows of ATP 
Synthase Dimers in Native Mitochondrial Inner Membranes. Biophys J 
2007;93:2870-2876. 

[10] Mastrangelo IA, Bezanilla M, Hansrna PK, Hough PV, Hansma HG. 
Structures of large T antigen at the origin of SV40 DNA replication by 
atomic force microscopy. Biophys J 1994;66:293-298. 

[11] Roos WH, Bruisma R, Wuite GJL. Physical virology. Nat Phys 
2010;6:733-743. 

[12] Cieplak M, Robbins MO. Nanoindentation of 35 virus capsids in a molec¬ 
ular model: Relating Mechanical Properties to Structure. PLOS ONE 
2013;8:e63640. 

[13] Lee EH, Gao M, Pinotsis N, Wilmanns M, Schulten K. Mechanical 
strength of the titin ZlZ2-telethonin complex. Structure 2006; 14:497509. 

[14] Bertz M, Wilmanns M, Rief M. The titin-telethonin complex is a di¬ 
rected, superstable molecular bond in the muscle Z-disk. Proc Natl Acad 
Sci USA 2009;106:13307-13310. 

[15] Dima RI, Joshi H. Probing the origin of tubulin rigidity with molecular 
simulations. Proc Natl Acad Sci USA 2008;105:15743-15748. 


22 


[16] Sikora M, Cieplak M. Mechanical stability of multidomain proteins and 
novel mechanical clamps. Proteins: Struct Funct Bioinf 2010;79:1786- 
1799. 

[17] Janowski R, Kozak M, Janowska E, Grzonka Z, Grubb A, Abraham- 
son M, Jaskolski M. Human cystatin C, an amyloidogenic protein, 
dimerizes through three-dimensional domain swapping. Nat Struct Biol 
2001 ;8:316320. 

[18] Zhmurov A, Brown AEX, Litvinov RI, Dima RI, Weisel JW, Barsegov 
V. Mechanism of Fibrin(ogen) forced unfolding. Structure 2011; 19:1615- 
1624. 

[19] Dietz H, Berkemeier F, Bertz M, Rief M. Anisotropic deformation re¬ 
sponse of single protein molecules. Proc Natl Acad Sci USA 2006; 
103:12724-12728. 

[20] Sikora M, Cieplak M. Formation of cystine slipknots in dimeric proteins. 
PLOS ONE 2013;8:e57443. 

[21] Sikora M, Cieplak M. Cystine plug and other novel mechanisms of large 
mechanical stability in dimeric proteins. Phys Rev Lett 2012;109:208101. 

[22] Cieplak M. Mechanostability of virus capsids and their proteins in 
structure-based models. In: Liwo A, editor. Computational methods to 
study the structure and dynamics of biomolecules and biomolecular pro¬ 
cesses - from bioinformatics to molecular quantum mechanics. Heidel¬ 
berg: Springer; 2014. p 295-315. 

[23] Velikovsky CA, Deng L, Chlewicki LK, Fernandez MM, Kumar V, Mar- 
iuzza RA. Structure of Natural Killer Receptor 2B4 Bound to CD48 


23 


Reveals Basis for Heterophilic Recognition in Signaling Lymphocyte Ac¬ 
tivation Molecule Family. Immunity 2007;27:572584. 

[24] Fuson KL, Ma L, Sutton RB, Oberhauser AF. The C2 Domains of Hu¬ 
man Synaptotagmin 1 Have Distinct Mechanical Properties. Biophys J 
2009;96:1083-1090. 

[25] Takahashi H, Shahin V, Henderson RM, Takeyasu K, Edwardson JM. 
Interaction of Synaptotagmin with Lipid Bilayers, Analyzed by Single- 
Molecule Force Spectroscopy. Biophys J 2010;99:2550-2558. 

[26] Cieplak M, Hoang TX. Universality classes in folding times of proteins. 
Biophys J 2003;84:475-488. 

[27] Sulkowska JI, Cieplak M. Mechanical stretching of proteins - a theoreti¬ 
cal survey of the Protein Data Bank. J Phys Cond Mat 2007; 19:283201. 

[28] Sulkowska JI, Cieplak M. Selection of optimal variants of Go-like models 
of proteins through studies of stretching. Biophys J 2008;95:3174-3191. 

[29] Sikora M, Sulkowska JI, Cieplak M. Mechanical strength of 17134 model 
proteins and cysteine slipknots. PLOS Comp Biol 2009;5:el000547. 

[30] Tsai J, Taylor R, Chothia C, Gerstein M. The packing density in pro¬ 
teins: Standard radii and volumes. J Mol Biol 1999;290:253-266. 

[31] Kwiecinska JI, Cieplak M. Chirality and protein folding. J Phys Cond 
Mat 2005;17:S1565-S1580. 

[32] Veitshans T, Klimov D, Thirumalai D. Protein folding kinet¬ 
ics: Timescales, pathways and energy landscapes in terms of sequence 
dependent properties. Folding and Design 1997;2:1-22. 


24 


[33] Szymczak P, Cieplak M. Stretching of proteins in a uniform how. J Chem 
Phys 2006; 125:164903. 

[34] Bell G. Models for the specihc adhesion of cells to cells. Science 
1978;200:618-627. 

[35] Evans E, Ritchie K. Dynamic Strength of Molecular Adhesion Bonds. 
Biophys J 1997;72:1541-1555. 

[36] Dndko OK, Hummer G, Szabo A. Intrinsic Rates and Activation Free 
Energies from Single-Molecule Pulling Experiments. Phys Rev Lett 
2006;96:108101. 

[37] Duan L, Zhmurov A, Barsegov V, Dima RI. Exploring the Mechanical 
Stability of the C2 Domains in Human Synaptotagmin 1. J Phys Chem 
B 2011;115:10133-10146. 

[38] Hu J, Lipowsky R, Weikl TR. Binding constants of membrane-anchored 
receptors and ligands depend strongly on the nanoscale roughness of 
membranes. Proc Natl Acad Sci USA 2013;110:15283-15288. 


25 


Tables and figures 


protocol and pathway 

P [e/A] 

q [e/A] 

ad [A] 

k 0 [1/ r] 

# trajectories 

N-N' pathway 1 

0.083 

2.66 

3.6 

2 • 10” 14 

373 

N-N' pathway 2 

0.089 

1.96 

3.4 

o 

1 

o 

92 

N-C' pathway 1 

0.089 

2.23 

3.4 

bO 

o 

1 

398 

N-C' pathway 2 

N/A 

N/A 

N/A 

N/A 

67 

C-N 7 

0.061 

1.74 

4.9 

9 • 10” 13 

420 

C-C' 

0.093 

2.52 

3.2 

2 • 10" 12 

420 

N-C pathway 1 

0.071 

2.20 

4.2 

O 

1 

231 

N-C pathway 2 

0.075 

2.43 

4.0 

2 • 10" 14 

239 


Table I: The fitting parameters p and q in Eq. (|3]l. and the resulting parame¬ 
ters ko and ad in Eq. (J4J , for the various pulling protocols and pathways. The 
differences in p and q between the protocols reflect the anisotropic response 
of the system to pulling. 
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2B4 



CD48 
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Figure 1: Left panel: structure of the CD48-2B4 complex. Proteins CD48 
(chain A) and 2B4 (chain B) are on the right (blue) and left (red), respec¬ 
tively. The residues that are involved in the CD48-2B4 binding are shown 
in the bond representation. The pulling simulations reveal three distinct 
patches of contacts at the CD48-2B4 interface, which are denoted by a (or¬ 
ange), b (purple), and c (green). The disulfide bridge between Cys 3 and 
Cys 100 in 2B4 is in yellow. The termini of CD48 and 2B4 are shown as 
spheres. We denote the CD48 termini by N and C, and the 2B4 termini by 
N' and C'. Right panel: the three patches of the inter-protein contacts are 
shown on CD48 on the side that forms the interface with 2B4. 
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Figure 2: Examples of simulation configurations of the CD48-2B4 complex 
in which shear (left panel) and tensile (right panel) mechanical clamps are 
formed. The arrows show the stretching directions, and the symbols (N and 
C in the left panel; N and N' in the right panel) lebel the termini to which the 
stretching forces are applied. The black dotted lines mark selected contacts 
forming the mechanical clamps. The residues involved are shown as spheres. 
The color code is the same as in Fig. [0 




0 40 80 120 0 40 80 120 

d [A] d [A] 

Figure 3: Examples of typical trajectories corresponding to two pathways 
for constant speed N-N' pulling. The top panels show the force-displacement 
curves. The bottom panels show the corresponding scenario diagrams in 
which the intra-protein contacts in CD48 are shown as black triangles whereas 
the contacts between CD48 and 2B4 are shown in orange, purple and green 
to match the three patches a, b and c defined in Fig. [TJ 
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0 50 100 150 0 50 100 150 

d [A] d [A] 

Figure 4: Similar to Fig. [3] but for pulling in the N-C' protocol. The arrows 
indicate the locations of the maximal force F max . 



0 18 36 54 0 12 24 36 

d [A] d [A] 

Figure 5: Similar to Fig. [3] but for the C-N' (left panels) and C-C' protocols 
(right panels). There are just single pathways in these protocols. 
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N-C 



Figure 6: Similar to Fig. [3] but for the N-C pulling protocol. The solid and 
dashed lines in the F-d curves correspond to the presence and absence of 
2B4, respectively. 
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d [A] d [A] 


Figure 7: Inter-residue distances r t j as a function of the spring displacement 
d. Left panel: Distances r l3 between residues i — 6 and j = 97, i = 7 and j = 
98, and i — 8 and j = 98, which represent three exemplary contacts between 
strands /3i and /3g in CD48. These residues are shown as blue spheres in the 
left-hand panel of Fig. [21 The moderate increase at d ~ 32 A corresponds 
to shear between f3\ and (3g. Right panel: Distances r VJ between residues 
i = 28 and j = 148, i = 35 and j = 203, and i = 88 and j = 145, which 
represent three selected contacts between CD48 and 2B4 in patches a, b and 
c, respectively. These contacts are shown in the right-hand panel of Fig. [21 
The steep increase at d ~ 125 A indicates action of the tensile mechanical 
clamps. 
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Figure 8: F max as a function of v v in the logarithmic scale for the indicated 
pulling protocols and pathways. The unit of speed is v 0 =l A/t. The dotted 
lines correspond to the relation F max = p log(v p /vo) + q with the fitting 
parameters p and q as specified in Table Q] 
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Figure 9: Pulling simulations of synaptotagmin 1. In the scenario diagram, 
the contacts within the separate C2A and C2B domains are shown as filled 
(blue) circles and empty (red) squares, respectively. The inter-domain con¬ 
tacts are represented by filled (cyan) diamonds. The highest peak in the F-d 
curve at d ~ 740 A is denoted by B. It is associated with rupturing the core 
of the C2B domain. The second highest peak at d ~ 60 A is denoted by A. 
It is associated with unraveling of fa and fa in the C2A domain. The minor 
force peak that is related to breaking the C2AB interface at d ~ 120 A is 
denoted by I. A typical simulation configuration at stage I is shown in the 
lower right corner. The dependences of F max on v p for the force peaks A and 
B are shown in the upper right panel (data averaged over a dozen trajectories 
at each of the pulling speeds). 
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Figure 10: Values of F max for the various protocols of pulling. The data above 
the diagonal correspond to v p = 5 x 10 -3 A/r. Those below the diagonal - 
to the pulling speed of 6 x 10” 6 A/r as obtained through extrapolation. A 
double entry indicates existence of two pathways. The entries just next to 
the diagonal are highlighted as they correspond to same-chain pulling. 
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Supporting Information 


This section contains four figures that support results reported in ” Unbinding 
and unfolding of adhesion protein complexes through stretching: interplay 
between shear and tensile mechanical clamps” by Bartosz Rozycki, Lukasz 
Mioduszewski and Marek Cieplak 
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Figure SI: The relative frequency of appearance of the first pathway, iq, as 
a function of the pulling speed v p for the N-N', N-C' and N-C protocols as 
indicated. 
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Figure S2: Separation distance d s as a function of the pulling speed v p in 
the logarithmic scale for the indicated pulling protocols and pathways. The 
dotted lines are shown to guide the eye. 
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Figure S3: The end-to-end distance at the time of CD48-2B4 separation, 
L s , as a function of the pulling speed v p in the logarithmic scale. The dotted 
lines are drawn to guide the eye. 
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Figure S4: Pulling simulations of synaptotagmin 1 with the inter-domain 
contact energies e' = 2 e (left) and e' = e/2 (right) at v p = 0.002 A/r. The top 
panels show typical F-d curves. The bottom panels show the corresponding 
scenario diagrams. Here, the contacts within the C2A domain are shown as 
filled (blue) circles, those within the C2B domain as empty (red) squares, 
and the inter-domain contacts as filled (cyan) diamonds. The two highest 
force peaks are indicated by dashed lines. 
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